Lattice Perturbation Theory by Langevin Dynamics 



Francesco Di Renzo 
Giuseppe Marchesini 
Paolo Marenzoni 
en ' Enrico Onofri 

Dipartimento di Fisica 
Universita di Parma, 
Viale delle Scienze 



< 



> 

o 
o 

oo 

o 

m 



X 



43100 Parma ( Italy 



Abstract 



We present here an application of the standard Langevin dynamics to the problem 
of perturbative expansions on the Lattice QCD. This method can be applied in the 
computation of the most general observables. In this work we will concentrate in 
particular on the computation of the perturbative terms of the lxl Wilson loop, up 
to fourth order. It is shown that a stochastic gauge fixing is a possible solution to the 
problem of divergent fluctuations which affect higher order coefficients. 



Since its introduction in 1981 by Parisi and Wu PJ, Langevin dynamics has been exten- 
sively used for Monte Carlo simulations. 

Basically it consists in a stochastic dynamical system on the field configuration space 
dictated by the general equation 

at d<p(x,t) +v[x ' th [i) 

where is the field, S[<p] the action and rj a Gaussian random noise satisfying to the nor- 
malization 

< 7](x, t)rf(x,t') >= 26 (x - x)6(t - t). (2) 



*We warmly thank P.Rossi, G. Parisi, G.C.Rossi, N.Cabibbo and C.Destri for interesting discussions and 
precious suggestions 



As a matter of fact, stochastic dynamics is devised in such a way that time averages on the 
noise converge to averages on the Gibbs measure 



< i J* dtO[<P(t)] >„- - I D0O[0]e- 5 ^. (3) 

To obtain from equation (1) a useful expression for computer simulations, one can take 
t discrete with time step dt = e: 

(j)(x,t n+ i) = <f>{x,t n ) - f x [<j>,T}], (4) 

where 

dS 

^ = e 7^ir~r\ + v^( z > *«) ( 5 ) 

0(f)(x,t n ) 

and now 77 is normalized by: 

< r](x h ti)T)(xj,tj) >= 25 XiXj 8 ti t r (6) 

In this discrete form, Langevin equation has to be regarded solely as an approximation of 
equation (1), valid only for e — ► 0. 

The method has been widely adopted as an alternative to Metropolis, Heat Bath and 
other Monte Carlo algorithms for scalar fields and gauge theories. In 1985 Batrouni et al. 
|| presented a new analysis for lattice fields theories. Given the standard Wilson action 

S=--j^ y ETr(U P + Ul), (7) 

where the sum is over the plaquettes P, in a four dimensional periodic lattice. Up are ordered 
products of the link gauge variables U^(x), which are SU(3) matrices. Here /x = 1, . . . 4 and 
x is a point of the lattice. Each configuration is then described by 4 x volume 3x3 complex 
matrices. 

For each link U = U^x), one adopts the evolution 

U(t n+1 )=e- F ^U(t n ), (8) 

where 



F <*■> - 1 



£ (c/ P - 4) - i £ zv(i7p-4) 



+ y/iH(t n ), (9) 



and if is a traceless antihermitian noise matrix with normalization given by 

< H ik (x,t)H lm (x ,t) > H = [5u5 km - -8ik8i m \8 xx '5 tt i. (10) 

Langevin approach was originally formulated for perturbation theory also on the contin- 
uum. What we present here is the application of this idea to compute the weak coupling 
expansion directly in the lattice. The problem is well known and has been considered by 



diagrammatic technique (see |4j]), which allows the calculation of the expansion coefficients 
up to g 4 (and in some cases g 6 ). Since gauge fields are written as 



U li = e* A -, 4 = -A„, TrAp — (11) 
where g is the coupling constant, the Langevin equation takes the following form: 

e 9< = e - F e gA ». (12) 
The fields can be expanded in series of g 

A, = Y.9 k Af. (13) 



In the same manner, recalling that (3 = 1/g 2 and imposing e = g 2 r, the drift (9) becomes; 



T 
12 



E E rr(?7p-4) 



+ gVfH(t n )=Y l 9 k F lk \ 



(14) 



The main point is to apply to equation (12) the Baker - Campbell - Hausdorff formula and 
to extract the contributions order by order in g. 

At present, we have implemented the simulation computing the evolutions for the gauge 
fields Ap up to fourth order in g. At this order, perturbative coefficients of many observables 
have been computed analytically. Thus, in order to check our lattice formulation of the 
above Langevin dynamics, the terms of the standard lxl plaquette have been measured 
(always to the order g 4 ), confirming the analytical results. 

While the original motivation of the Langevin approach was to make it possible to cal- 
culate in perturbation theory without fixing a gauge, it is known that some divergent fluc- 
tuations (averaging to zero) may plague high order terms. We confirm this phenomenon in 
the case of the plaquette expansion. We observe indeed that, from the third order in g, the 
errors associated to our observables grow in time, even if the mean value remains always 
stable around its known value. In higher orders, this spurious fluctuation may completely 
hide the signal. As a way out we have applied a technique which goes back to Zwanziger 
and has been more recently implemented on the lattice |J in the form of "stochastic gauge 
fixing" . We may think to a new source in the Langevin equation, responsible of a stochastic 
gauge fixing, so that the actual algorithm implemented is (for the gauge links) 



U' N (x) = e F W™W(x) 



U^ +1 (x) = e w ^ N ' x ^U' N (x)e~ w ^ N ' x ' IJ 'y 



(15) 



with 



w[U„,x) =«E, (A_„ [U^x) - Ul(x - fi)]) 
A^U v (x) = U„(x) - U v (x - n) 



traceless 



(16) 



a being a free parameter. 

As a matter of fact, the result is very impressive. We report, for example, the term in g A of 
the plaquette, measured (in figure 1) without gauge fixing, and (in figure 2) with the above 
gauge fixing. By means of this essential reduction of the noise, the result we obtained for 
the reported term, with r and a set to 0.02, is c 4 = 1.211 ± 0.003, through an average over 
2048 iterations. 

All the numerical experiments have been done on Connection Machine CM2, using CM 
Fortran, with a great contribution of the CMIS assembler. This kind of computation is 
indeed very expensive: 2.455 Gflops are needed to complete a single Langevin iteration on 
a lattice size of 8 4 , at the fourth order, with the measure of the corresponding coefficient of 
the plaquette. For the same lattice size, the program uses about 10 MBytes of memory. 

The work is in progress to go to higher orders and test a wide class of observables of 
interest for Lattice gauge Theories. 
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